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Under which conditions do the electrical transport properties of one-dimensional (ID) carbon 
nanotubes (CNTs) and 2D graphene become equivalent? We have performed atomistic calculations 
of the phonon-limited electrical mobility in graphene and in a wide range of CNTs of different types 
to address this issue. The theoretical study is based on a tight-binding method and a force-constant 
model from which all possible electron-phonon couplings are computed. The electrical resistivity of 
graphene is found in very good agreement with experiments performed at high carrier density. A 
common methodology is applied to study the transition from ID to 2D by considering CNTs with 
diameter up to 16 nm. It is found that the mobility in CNTs of increasing diameter converges to 
the same value, the mobility in graphene. This convergence is much faster at high temperature and 
high carrier density. For small-diameter CNTs, the mobility strongly depends on chirality, diameter, 
and existence of a bandgap. 


I. INTRODUCTION 

Graphene and carbon nanotubes (CNTs) are respec¬ 
tively two-dimensional (2D) and one-dimensional (ID) 
allotropes of pure sp 2 carbon. Their remarkable physical 
and chemical properties have been studied extensively^— 
but investigations of the electrical transport in graphene 
and CNTs are usually done separately. Of course, it 
is known that graphene is a semi-metal characterized 
by Dirac cones at the charge neutrality point while 
CNTs can either be metallic or semiconducting depend¬ 
ing on their chirality. Carrier mobilities in graphene 
and CNTs may reach very high values, which make 
them promising for high-frequency applications. How¬ 
ever, record values differ considerably between graphene 
(> 10 5 cm 2 /V/s )££ and semiconducting CNTs (> 
10 3 cm 2 /V/s)— r— Furthermore, one should in principle 
recover the transport properties of graphene by increas¬ 
ing the diameter of metallic or semiconducting CNTs, at 
least in situations where this transport is only limited by 
intrinsic processes like phonon scattering. Quite supris- 
ingly, this transition of the transport properties from ID 
to 2D has never been carefully investigated even if there 
are numerous studies on the electron-phonon coupling in 
the two materials — — 

In this work, we study the evolution of the phonon- 
limited mobility of carriers from CNTs to graphene, i.e., 
from ID to 2D. We calculate the electrical mobility (or 
resistivity) using a common methodology, enabling a 
direct comparison between graphene and semiconduct¬ 
ing/metallic CNTs of varying diameters. We focus on 
phonon scattering not only because it is a limiting mech¬ 
anism at high carrier density and high temperature^ 
but also because it is intrinsic to the materials, at vari¬ 
ance with scattering by impurities or by surface optical 
phonons of nearby oxides, which are dependent on the ge¬ 
ometry, on the quality of the materials and on the nature 


of the environment. We also leave aside carbon nanorib¬ 
bons in this study, since the possible existence of edge 
states does not allow for a straightforward analysis of 
the ID to 2D transistion. 

In the following, we present fully atomistic calculations 
of the phonon-limited mobility in graphene and CNTs 
with diameter up to 16 nm. We show that the mobil¬ 
ity in CNTs tends to the limit in graphene at large di¬ 
ameter (sometimes above 10 nm), but that the route to 
it strongly depends on size, chirality, and temperature. 
These behaviors are explained by specific features of the 
band structure and of the electron-phonon coupling. 

II. METHODOLOGY 
A. Justification of the methodology 

A strong motivation for this study is that it sheds light 
on the effects of confinement and dimensionality upon the 
electron/phonon band structures, the electron-phonon 
coupling and the transport. For example, the scattering 
of carriers in ID semiconductor nanostructures such as 
nanowires is merely limited to back-scattering when only 
one electronic band is populated. This tends to reduce 
the scattering probability. At the same time, confinement 
leads to an enhancement of the electron-phonon coupling, 
which is stronger in ID than in 2D (and 3D).— 7 In this 
context, it is particularly interesting to investigate these 
effects in ID and 2D carbon allotropes since their band 
structures considerably differ from those of usual semi¬ 
conductor nanostructures. In addition, it is important to 
understand how (and when) the phonon-limited carrier 
mobility in CNTs of increasing diameter approaches the 
limit in graphene as a function of temperature, carrier 
density, chirality, and metallic or semiconducting charac¬ 
ter. 
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There are already numerous theoretical studies on the 
phonon-limited transport in graphen e 15 ' 17-19 ' 22 ' 23 ' 28 and 
CNT s 2 ' 11 ' 29- — but it is only in recent works on graphene 
that the couplings to all phonons have been calculated 
from first-principles without extracting parameters from 
experiments . 20 ' 22 ' 23 From these works, we learn that it is 
essential to consider many different phonon modes like 
longitudinal acoustic (LA), transverse acoustic (TA) and 
optical phonons to predict the temperature- and density- 
dependent resistivity of graphene. We deduce that it is 
likely the case for CNTs, especially because the k selec¬ 
tion rule is broken in the direction perpendicular to the 
CNT axis, in particular in small-diameter CNTs. 

Therefore it is essential to consider all possible 
electron-phonon scattering processes to predict the 
phonon-limited mobility in CNTs. In order to catch the 
1D-2D transition, it is also necessary to go well beyond 
previous calculations that were limited to small-diameter 
tubes. However, first-principles calculations cannot be 
performed on large-diameter CNTs because the compu¬ 
tational cost increases dramatically with the number of 
atoms per unit cell. In this work, we present calcula¬ 
tions combining tight-binding (TB) and force-constant 
models for electrons and phonons, respectively. The pa¬ 
rameters of these models were refined on the latest first- 
principles calculations; 23 ' 36 The low-field mobility is com¬ 
puted taking all electron and phonon bands into account, 
the carriers being coupled to all possible phonons, in¬ 
cluding intra- and inter-subband scattering. We show 
that this approach gives results for graphene in excellent 
agreement with first-principles calculation s 23 ’ 36 and with 
experiments ; 24 ' 37 The same methodology is then applied 
to CNTs. In the following, we present results for p-type 
graphene and CNTs, i.e., for the transport of holes. Very 
similar results are found for electrons. For the sake of 
comparison, carrier densities are given as equivalent sur¬ 
face densities (in cm -2 ). For a CNT of diameter d , the ID 
and 2D carrier densities are related by 712 d = nm/(jrd). 


B. Tight-binding Hamiltonian 

The electronic states of graphene and CNTs are writ¬ 
ten in the basis of the p z orbitals where z is the axis per¬ 
pendicular to the lattice. Interactions are restricted to 
first-nearest-neighbors (INNs). The Hamiltonian matrix 
is therefore defined by two quantities, the onsite energy 
E p and the INN hopping term V PP n that both depend on 
the atomic displacements induced by the phonons. The 
variations of V pp „ are governed by a power law^ 

V =V° (*>Y~ (1) 

VppiT — v pp7r I ^ I J 

in which d and do are the actual and equilibrium bond 
lengths between INN atoms, respectively. The param¬ 
eters in Eq. JTJ were set based on the GW calculation 
data in Ref. [36J, i.e. do = 1.42 A, V p = 3.0 eV, (which 


gives 10 6 m/s for the Fermi velocity,) and n pp7r = 3.47, 
(which fits the canonical acoustic gauge field, 5.52 eV). 

The variations of the onsite term E p with atomic 
displacements were ignored in most previous TB 
studie s 15 ' 19 ’ 28 ' 29 arguing that they give little contribu¬ 
tion to the deformation potentials. In our case, we follow 
the arguments of Ref. [39j and introduce in E p a term 
proportional to the average of the relative bond length 
variation 



p 


E 



( 2 ) 


where the sum is over the three INN atoms. The Dirac 
point at d = do is taken as the origin of the elec¬ 
tronic energies. The parameter a p in Eq. [2] is set to 
—4.162 eV, and has been extracted from Density Func¬ 
tional Theory (DFT) calculations of the band structure 
of graphene for different cell sizes. These DFT calcu¬ 
lations were performed with the BigDFT— code, us¬ 
ing Norm-Conserving Pseudopotentials^ and the Gen¬ 
eralized Gradient Approximation^ for the exchange and 
correlation potential. BigDFT makes use of a system¬ 
atic real-space wavelet basis, and can handle surface-like 
boundary conditions^ which allows to define a common 
energy reference for all cell sizes and enables a direct 
comparison of energy levels; 44 

The changes in V/ p7r due to bond length variations 
(Eq. [J) contribute to the coupling of electrons to both 
LA and TA phonons, whereas those in E p (Eq. [2]) are 
only relevant for LA phonons. Their respective influ¬ 
ence on the mobility was heavily debate d 10 ' 11 ' 15 ' 19 ' 28 ' 29 
but it is now clear that both of them must be taken into 
account.— 

In the case of CNTs, the hopping terms between INN 
p z orbitals slightly vary from one bond to another due to 
the curvature . 2 ' 45 ' 46 As a consequence, there is a small 
curvature-induced gap in ’’metallic” zigzag CNTs scaling 
as 1/d 2 . For example, we obtain gaps of 80, 45, and 29 
meV in (9,0), (12,0), and (15,0) CNTs, respectively, in 
excellent agreement with experiments.— For symmetry 
reasons, armchair CNTs always preserve their metallic 
character . 2 ' 45 ' 46 The a — n hybridization induced by the 
curvature^ is neglected in the present calculations. 


C. Phonons 

For the calculation of the phonons, the dynamical ma¬ 
trix is built by a fourth-nearest-neighbour (4NN) force- 
constant model as originally fitted to experimental data 
by Jishi et al.— and later fitted to ab-initio calcula¬ 
tions and enhanced by off-diagonal force-constants for 
the 2NN.— For the current calculations, we refitted the 
model, taking also the off-diagonal elements of the 4NN 
into account. The importance of the 4NN off-diagonal el¬ 
ements had already been noted before in ab-initio phonon 
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calculations— While the phonon-dispersion itself can be 
fitted very well without the off-diagonal terms, they have 
an important impact on the phonon modes: they are nec¬ 
essary to obtain the correct amplitude for the admixture 
of optical components to the acoustic phonon modes at 
non-vanishing wave-vector. This admixture considerably 
influences the contribution of the acoustic modes to the 
mobility— For example, the ratio of the acoustic phonon 
deformation potential with and without optical phonon 
mixing is about 64.2% in our model, which is compara¬ 
ble to the one in Ref. [36|. The parameters of our force- 
constant model and the dispersion obtained from it are 
given in Appendix lAl 


The electron-phonon scattering rates are calculated for 
all electronic bands within at least 5fcsT+200 meV of the 
Fermi level, and for all phonon modes. 

E. Mobility 

The low field mobility p is obtained by the resolution 
of the Boltzmann transport equation in the stationary 
regime. Under the application of a constant electric Held 
F, the distribution function in the state |k, b) is given to 
the first-order in F by f b (k) = /°(£k,b) + eFgb(k) where 
f° is the Fermi-Dirac distribution function. In CNTs, 
< 7 &(k) is solution of the following equations: 


D. Electron-phonon coupling 

The same parameters are used for the calculations of 
the electron and phonon band structures of CNTs and 
graphene - there is no additional parameter. 

Following Ref. [2^, the matrix element for the 

transition of an electron from an initial state |k, b) to a 
final state |k', b') (k' = k + q) after emission of a phonon 
I — q, j) (same formula for absorption) is given by 


/ <7b(k){ Wkb,k'6'[l — /°(-Ek',&')] + Wk'&',kb/°(-E ; k',&')} 

v 


— 9b'{k'){Wk'b',kb 
'df_ 
dE 


2n ziu 
= ~ Vb W 


[1 - f°(E ktb )] + W kb ^f°(E k , b )}dk' 

) ( 5 ) 


where v&(k) = H~ 1 dE k _ b /dk is the group velocity along 
the electric field and L is the length of the CNT. The 
mobility is then given by 
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where k and q are the wave vectors of the electron and 
phonon, respectively, while b and j are the indexes of the 
electron band and phonon mode. C^' b is an eigenvec¬ 
tor element of the electronic Hamiltonian H at equilib¬ 
rium, where f) (or a) is the atom index in the unit cell. 
cf>(r — R m|S ) is the p z atomic orbital centered on atom /? 

in the unit cell to. (q) is an eigenvector element of 
phonon state |q, j) and Wy (q) is the corresponding eigen- 
frequency. Rod is the i component (a ’,y,z) of vector R 0a , 
N is the number of Wigner-Seitz unit cells, and M is the 
mass of a carbon atom. The transition rate is given by 
Fermi’s golden rule, 


_ I , , 2 

H'kh.k'b' = |-^k,6 b 

0 

X {uq ) jA(F/k / ,6 / E kb Hujj( q)) 

+ [^q ,j + l]^(£ , k',b' — E kb + hu!j( q))} (4) 


where E k b is the energy of |k, b) and n qj - is the equi¬ 
librium phonon occupation number (Bose-Einstein dis¬ 
tribution) . 


_ _ b J g b {k)v b {k)dk 

M E b ff°(E Kb )dk ■ 1 j 

Equations ([5][6]) remain valid for graphene using the 
substitutions dk' —> d 2 k' and 27 t/L —> 4:Tt 2 /S where S is 
the sample surface. 

The linearized Boltzmann transport equation is solved 
exactly. Brillouin zone integrations are performed on a 
non-homogeneous ID grid for CNTs and on an triangular 
mesh for graphene (e.g., 8600 fc-points in ID and 17124 
triangles in 2D for a surface carrier density of 13.6 x 
10 12 cm” 2 and T = 300 K). 

III. RESULTS AND DISCUSSION 
A. Application to graphene 

Figure [T] presents the electrical resistivity p of graphene 
calculated versus temperature T, for several carrier 
densities U 2 d- In accordance with early theoretical 
predictions^ and experiments,— the resistivity is pro¬ 
portional to T 4 at low temperature, while at higher tem¬ 
perature it varies linearly with T (the threshold strongly 
depends on U 2 d)- In the T 4 regime, short wavelength 
acoustic phonons are frozen out, restricting scattering 
processes to small scattering angles—*? 4 In agreement 
with DFT calculations,— both LA and TA phonons con¬ 
tribute to the scattering in the T 4 and T regimes. We also 
confirm that: i) in the T regime and above, the resistiv¬ 
ity does not depend on n 2 flr’ 24 ii) at room temperature 
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FIG. 1. Comparison of the resistivity of graphene versus tem¬ 
perature for different hole densities, (a) Present work (solid 
lines) and experimental data (dotted lines), which correspond 
to the temperature dependent part of the resistivity mea¬ 
sured in Ref. l2ll . (b) Present work (solid lines) and DFT 
calculations.— 


and above, the resistivity is enhanced by optical phonon 
scattering 20 ^ so that it deviates from the T trend. 

In Ref. H3, the resistivity of graphene was measured at 
high carrier density (n 2 D > 10 13 cm -3 ). Its temperature 
dependent component is displayed in Fig. |T] and is com¬ 
pared to our predictions. The agreement between theory 
and experiments is fairly good. As shown in Appendix iBl 
in this range of carrier density, the resistivity is mainly 
limited by (intrinsic) phonon scattering, justifying the 
comparison. Also, the agreement is excellent with the 
DFT calculations of Ref. O (Fig. [TJd). 

In Appendix [Bj we show that our predicted resistivi¬ 
ties also agree with the experimental data of Ref. H3 mea¬ 
sured at lower carrier density. In that case, we have to 
include scattering by the surface optical (SO) phonons of 
the substrate. We also demonstrate in Appendix IBl that 
scattering by intrinsic phonons becomes more efficient 
than scattering by SO phonons not only at low temper¬ 
ature because it involves lower energy phonons, but also 
at high temperature and high carrier density where SO 
phonon potentials are efficiently screened. 


T = 300K, ri = 1 x 10 12 cm -2 





FIG. 2. Plionon-limited mobility in CNTs with different chi¬ 
rality (armchair (n, n), zigzag (n, 0), chiral (2 n,n) and chi¬ 
ral (4n, n)) compared to graphene (green line) versus (a) di¬ 
ameter at 300 K and ri 2 D = 10 12 cm -2 , (b) carrier den¬ 
sity at 300 K and d « 3.5 nm and, (c) temperature at 
ri 2 d = 13.6 x 10 12 cm -2 and d « 3.5 nm. 


B. Carrier mobility in CNTs 


We consider now zigzag (n,0), armchair ( n,n ), chiral 
(2n, n ) and chiral (4n, n ) CNTs. (n, n ) and (4n, n ) CNTs 
are always metallic. Zigzag (n, 0) and chiral (2n, n) CNTs 
are metallic when n is a multiple of 3, semiconducting 
otherwise. As discussed above, a gap is induced by the 
curvature in small-diameter metallic (n, 0) CNTs. The 
electronic band structures of CNTs (Appendix [Cl) are in¬ 
timately related to that of graphene*^ In the 2D k space 
of graphene, the ID k vectors allowed for CNTs form a 
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bundle of parallel lines along the tube direction, the sep¬ 
aration between the lines being inversely proportional to 
the CNT diameter. This results in CNT band struc¬ 
tures composed of successive subbands, one for each line. 
In these conditions, since carrier transport takes place 
within an energy window of a few ksT around the Fermi 
surface, we anticipate that the transport properties of 
CNTs can reach those of graphene in three ways, namely 
by increasing 1) the tube diameter, 2) the carrier density, 
or 3) the temperature. In the three cases, the number of 
subbands in the transport energy window will increase 
in such a way that the ID transport will progressively 
turn into 2D transport - the quantification of k perpen¬ 
dicular to the tube axis becoming irrelevant. Figure [2] 
illustrates the evolution of the mobility along the three 
ways. It shows that the mobility in CNTs actually tends 
to its limit in graphene but with very different behaviors 
depending on the nature of the CNTs. The reasons why 
are discussed below. 


1. Mobility versus diameter 

The carrier mobility in CNTs calculated for ri 2 D = 
10 12 cm -2 and T = 300 K is plotted as a function of 
the diameter d in Fig. [2ji. The mobility is consider¬ 
ably smaller in semiconducting CNTs with d < 10 nm 
than in metallic ones, and is also much smaller than in 
graphene for the same carrier density. This is mainly 
due to the non-zero mass of the carriers since the bands 
are parabolic near their extrema. For d < 4 nm, the 
mobility in zigzag semiconducting CNTs scales with di¬ 
ameter as d 21 , in agreement with early predictions^ 
and measurements.- 7 An approximately quadratic rela¬ 
tionship is also obtained for small-diameter (2n, n) semi¬ 
conducting CNTs, as reported for different chiralities^ 
The important reduction of the mobility in semiconduct¬ 
ing CNTs with decreasing size is due to the increase of 
the effective mass^ and to the enhancement of electron- 
phonon coupling, in particular with radial breathing 
modes.- ’ 29 ’ 31 ’ 32 The strong enhancement of the electron- 
phonon coupling in small-diameter ID conductors seems 
to be a general trend, as shown for example for electrons 
in thin Si nanowires ^ - i— For d > 10 nm, the mobility 
in semiconducting zigzag CNTs saturates at a constant 
value, however smaller than in graphene for these par¬ 
ticular T and ri 2 d- Quite generally, we have found that 
plateaus appear in mobility curves when the Fermi level 
lies between two subbands and the upper subband does 
not contribute to transport. 

For metallic CNTs, the behavior of the mobility versus 
diameter is very different (Fig. [2ji), as expected from the 
presence of linearly dispersive bands (Appendix [Cl). 31 ’ 52 
At large diameter, the mobility is close to the value 
in graphene as several bands are populated, including 
parabolic bands that contribute to the diffusive trans¬ 
port. Interestingly, for U 2 d = 10 12 cm -2 and T = 300 K, 
the mobility in metallic zigzag CNTs saturates at the 


same plateau as in semiconducting ones for 12 < d < 
16 nm. 

At smaller diameter (4 < d < 12 nm), only linearly 
dispersive bands are occupied in metallic CNTs (Ap¬ 
pendix [Cj, and the mobility increases because the num¬ 
ber of allowed scattering processes decreases. The mo¬ 
bility is larger in armchair than in zigzag CNTs because 
there is no coupling to the LA phonons in the linear 
bands of armchair CNTs for symmetry reasons^. The 
total scattering by acoustic phonons is, therefore, much 
smaller in armchair CNTs (see Sect. IIII B 41 below!. The 
difference between metallic armchair and zigzag CNTs 
is progressively reduced when the temperature or car¬ 
rier density are increased (Fig. [2}j,c) , since carriers start 
to occupy parabolic bands in which the scattering by LA 
phonons is allowed, in both armchair and zigzag CNTs^i 

In metallic zigzag CNTs with small diameter (d < 
3 nm), the mobility suddenly decreases because the 
electron-phonon interaction considerably increases, in 
particular with LA and radial breathing modest as in 
semiconducting CNTs. For very small diameters, the 
bands near the neutrality point become parabolic due to 
the curvature-induced gap but the effect on the mobility 
is small for all carrier densities considered here. 

The mobility in chiral metallic CNTs [(2n, n) or 
(4n, n)\ follows the same trends as in zigzag metallic 
CNTs versus diameter, carrier density or temperature. 
Similarly, the (2 n,n) chiral CNTs show the same be¬ 
havior as zigzag semiconducting CNTs. Therefore, chi¬ 
ral and achiral CNTs of the same nature (semiconduct¬ 
ing/metallic) exhibit similar phonon-limited transport 

properties^ 


2. Mobility versus carrier density 

We discuss now the evolution of the carrier mobility 
as a function of n 2 Dj the other parameters (T = 300 K, 
d « 3.5 nm) being kept constant (Fig. [2Jj). For car¬ 
rier densities of the order of 10 12 cm~ 2 , the mobilities of 
the different types of CNTs span more than one order of 
magnitude. Two trends can be highlighted when ri 2 D is 
increased. First, the mobility decreases, almost exactly 
like l/ri 2 D as in the case of graphene (see the plot of 
/r x n- 2 d versus ri 2 D in Appendix [D]). The well-known in¬ 
dependence of /i x n, 2 D (or the resistivity) on the Fermi 
energy or equivalently on the carrier density in graphene 
at 300 is approximately recovered in CNTs at high 

carrier density. Second, above ~ 2 x 10 13 cm -2 , the mo¬ 
bilities of CNTs and graphene tend to converge to the 
same value. In that case, the integration over many sub¬ 
bands smooths out the effects resulting from the ID char¬ 
acter of the CNTs. For n 2 D of the order of 10 12 cm -2 , 
the mobility is higher in metallic CNTs than in graphene 
but the situation is reversed when ri 2 D reaches a certain 
threshold (for example ~ 7 x 10 12 cm -2 for armchair 
CNTs). In fact, the mobility becomes smaller in metal¬ 
lic CNTs than in graphene when the parabolic subbands 
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start to be occupied (Appendix [C]). The opposite situa¬ 
tion arises when the Fermi level only crosses the linear 
bands. 


3. Mobility versus temperature 

The variations of the mobility with temperature at 
fixed carrier density (10 12 cm~ 2 ) and diameter (3.5 nm) 
are plotted in Fig. [5J;. At high temperature, the mobility 
in CNTs decreases roughly like 1 /T@ where the exponent 
P is close to unity like in graphene but slightly differs from 
case to case. However, for this particular diameter and 
carrier density, the spread of the calculated mobilities is 
still important at 500 K. This shows that the broaden¬ 
ing and the averaging induced by the temperature are 
not sufficient to smooth out ID band structure effects 
at this diameter. As already discussed for graphene, the 
exponent /3 differs from unity because of the increasing 
role of optical phonons at high temperature. However, 
we show in the next section that the relative importance 
of optical phonons at high T strongly varies from one 
CNT to another, being huge in armchair metallic CNTs 
and remaining modest in semiconducting CNTs [zigzag 
or (2 n, n)]. 

4- Role of the different phonons 

In this section, we clarify the relative importance of the 
different phonons for each type of CNTs, and we compare 
with the case of graphene. 

Figure [3] shows the ratio between the mobility cal¬ 
culated by considering phonons with energy up to hu 
and the mobility including all phonons. In graphene, 
this ratio has a plateau around 1.5 from 30 to 130 meV 
(green line on Fig. [3))). This is consistent with previ¬ 
ous studies^ which reported that high-energy phonons 
contribute to 30% of the resistivity in graphene at room 
temperature. Similar plateaus are found in CNTs but 
at different values: 2.5 for armchair CNTs, 1.25 for 
other metallic CNTs, and 1.05 for semiconducting CNTs 
(d ss 3.5 nm). This means that the relative impact of high 
energy phonons strongly depends on the type of CNT. 
In armchair CNTs, this impact is huge - about 60% of 
the resistivity - because the coupling to LA phonons is 
forbidden in the linear bands and therefore the relative 
contribution of the optical phonons is more important. 
The latter is about 20% in other metallic CNTs, and 
about 5% in semiconducting CNTs, since the coupling to 
acoustic phonons is typically strong in parabolic bands. 

Figured] shows that 33% of the electrical resistivity of 
graphene at 300 K comes from the scattering by high- 
energy (>130 meV) phonons, in agreement with recent 
ab-initio calculations.— This ratio reaches 60% in 3.5 nm 
armchair metallic CNTs where the coupling to the LA 
phonons is weak. On the contrary, this ratio is only 5% in 
3.5 nm semiconducting CNTs because the scattering by 





FIG. 3. Ratio between the mobility calculated by consider¬ 
ing phonons with energy up to hw and the mobility calcu¬ 
lated with all phonons (a) for armchair CNTs; (b) for zigzag, 
(2n, n), (4n, n) metallic CNTs and graphene; (c) for zigzag, 
(2n, n) semiconducting CNTs. d « 3.5 nm, n = 10 12 cm -2 . 



FIG. 4. Ratio between the mobility calculated by considering 
phonons with energy below 130 meV and the mobility cal¬ 
culated with all phonons versus temperature, d « 3.5 nm, 
n = 10 12 cm' 2 . 
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FIG. 5. Phonon-limited mobility in CNTs versus diameter 
for T = 100 K (dashed lines), 300 K (dash-dotted lines) and 
500 K (dotted lines). Carrier density: (a) 10 12 cm -2 , (b) 
13.6 x 10 12 cm” 2 . The horizontal lines indicate the mobility 
for graphene. The horizontal scale is smaller in (b) than in 
(a), highlighting the faster convergence of all the mobilities 
at high carrier density. 


low-energy acoustic phonons is very strong and the mo¬ 
bility is low. The coupling to optical phonons also plays 
an important role in high-field transport in graphene and 
CNTs. 53 ' 54 The impact of high-energy phonons increases 
with the temperature (Fig. [4]). The contribution of these 
phonons to the resistivity reaches 10% at about 180, 210, 
250 and 350 K in armchair CNTs, graphene, metallic 
CNTs and semiconducting CNTs, respectively. 


5. Mobility from CNTs to graphene 

Figure [5] illustrates differently how the mobility in 
CNTs converges towards the value in graphene when the 
diameter is increased. The convergence is quite fast at 
high carrier density (13.6 x 10 12 cm -2 ) and high tem¬ 
perature (T > 300 K, Fig. [5}j). For example, at 500 K, 
the mobility is basically the same in the three kinds of 
investigated CNTs as in graphene for d ~ 5 nm. The 
transition from ID to 2D transport takes place above a 


certain diameter which strongly depends on temperature 
and carrier density. Below this threshold, the transport 
properties strongly depend on the chirality through its 
effect on the band structure. 


IV. PERSPECTIVES 

Experimentally, the diameter of single-walled CNTs 
which are usually grown is typically between 0.7 and 
4 nm. However, a recent work reported the synthesis of 
single-walled CNTs in the 5-10 nm diameter range (5% of 
the CNTs even have a diameter above 10 nm).— There¬ 
fore we believe that the experimental observation of the 
1D-2D transition is possible, in particular using the same 
approach allowing to reach ultrahigh carrier densities up 
to 4 x 10 14 cm -2 in graphene— 

Another interesting perspective of the present work 
would be to consider the effects of electron-electron inter¬ 
actions on the 1D-2D transition in CNTs. As discussed 
above, the transport properties in graphene at high car¬ 
rier density and temperature (> 20 K) are unambigu¬ 
ously dominated by clectron-phonon scattering. How¬ 
ever, theoretica l 56 ' 57 and experimental— studies show 
that small-diameter metallic CNTs exhibit Luttinger- 
liquid behaviour characterized by low-energy collective 
excitations of the electrons. In that case, the tunneling 
of carriers at energies near the Fermi level is strongly 
suppressed, so that the conductance increases as a power 
law with respect to temperature^ On the contrary, the 
conductance in semiconducting CNTs strongly decreases 
with temperature, as expected for diffusive transport 
limited by phonons J All these works deal with small- 
diameter CNTs at low carrier density. Therefore it would 
be extremely interesting to extend experimental stud¬ 
ies on CNTs to larger diameter and higher carrier den¬ 
sity. The electron-phonon interactions should indeed 
overcome electron-electron interactions when the num¬ 
ber of populated subbands increases, in particular for 
parabolic bands. In these conditions, the carrier mobil¬ 
ity in metallic CNTs should follow our predictions above 
certain thresholds that must be found. 


V. CONCLUSION 

It is shown that the phonon-limited mobility in ID 
CNTs approaches that of 2D graphene continuously by 
increasing the size of CNTs, the carrier density, or the 
temperature. The physics of this transition has been 
studied using atomistic calculations combining a tight- 
binding model for electrons, a force-constant model for 
phonons, the computation of all the electron-phonon cou¬ 
plings, and a full resolution of the Boltzmann trans¬ 
port equation. This approach gives carrier mobility in 
graphene in excellent agreement with experiments^ and 
DFT calculations^ and therefore can be used to pre¬ 
dict the mobility in CNTs in a wide range of configura- 
























tions. The mobility in CNTs can be higher or smaller 
than in graphene depending on chirality, diameter, car¬ 
rier density or temperature but converges to the same 
value above varying thresholds. This ID to 2D transi¬ 
tion takes place when the number of subbands situated 
in the transport energy window is sufficiently large to 
smooth out the effect of the ID confinement on the band 
structure and on the electron-phonon coupling. 

Appendix A: Parameters for the 4NN force constant 

model 

We have refitted the force constant method of Ref. [5C] 
to an ab-initio calculation of the phonon-dispersion of 
graphene using the ABINIT code4^ The electronic struc¬ 
ture is calculated in the local-density approximation 
(LDA) using a regular 60x60 reciprocal mesh in the first 
Brillouin zone and an energy cut-off of 35 Ha using a LDA 
functional. A thermal Fermi-Dirac smearing of 0.002 Ha 
is employed. We find the optimized cell parameter to 
be 4.631 A. The dynamical matrices were calculated us¬ 
ing density functional perturbation theory (DFPT) on a 
30x30 q-mesh. Since LDA overbinds, i.e., phonon fre¬ 
quencies have the tendency to be slightly too high, a 
scaling factor is used such that the phonon frequencies 
of the LO/TO mode at Gamma match the experimental 

valued 

The (real-space) force constants between two particu¬ 
lar atoms a and b are defined as the second derivatives 
of the total energy of the system with respect to the dis¬ 
placements of atom a in direction i and of atom b in 
direction j: 

a = Ai£_ (ah 

In local coordinates, direction aq is along the line con¬ 
necting the two atoms, X 2 is perpendicular to this line in 
the plane of the graphene sheet, £3 is the out-of-plane 
direction. In this local reference frame we define the 
longitudinal forces, (< j) l n ), transverse in-plane (</>((), and 
transverse-out-of-plane (</>(/) forces that act on a partic¬ 
ular atom when its nth nearest-neighbor is displaced. In 
the conventional 4NN-force constant model, only these 
“diagonal” terms are fitted. In our model, we include 
also the “off-diagonal” coupling between the longitudinal 
direction and the transverse-in-plane direction ( e\\( U and 
&n/ 1 ). The force-constant matrix for the interaction be¬ 
tween two atoms in the local reference frame thus reads: 

(€ ^ 0 \ 

c; = et /l 0 (A2) 

V 0 0 4L°J 

The off-diagonal force constants eK U and Cn' 1 obey the 
following relations 



FIG. 6. Phonon dispersion of graphene. The red lines are the 
frequencies calculated using the 4NN force constant model, 
and grey lines are DFT-LDA calculations using ABINIT. 


l/ti _ l/ti _ p. l/ti _ ti/l 

un = u/i = 0 e i/u = un (A3 ) 

During the fitting, we have noticed that the off- 
diagonal terms of the second and the fourth nearest 
neighbor (1 e l / tl and e l / u ) are essential. While a fairly 
decent fit of the phonon-dispersion alone can be achieved 
without these off-diagonal terms, the correct ratio of the 
amplitudes of optical and acoustic phonon components 
in the transverse and longitudinal acoustic branches at 
q 0 can only be achieved with the inclusion of the off- 
diagonal forces. By using the parameters in Table |U 
the ab-initio ratio was reproduced to a very good degree. 
The comparison of the phonon dispersion relations be¬ 
tween the ab-initio calculation and the 4NN force con¬ 
stant model is shown in Fig. [ 6 j We note that the 
agreement is good but the 4NN model cannot reproduce 
the Kohn anomalies in the two highest-optical branches. 
This would require the inclusion of many more distant 
neighbor interactions in the model (even infinitely many, 
if one wants to reproduce the kink^i). However, this does 
not seem to be necessary for our present study since the 
4NN model gives electron-phonon scattering rates and 
carrier mobilities in excellent agreement with ab-initio 
calculations for graphene^ 


Appendix B: Resistivity of graphene including 
surface optical phonon scattering 

The resistivity of graphene was measured in Ref. [13 for 
carrier densities varying from 1 to 3 x 10 12 cm -2 . In the 
following, we show that this resistivity consists of three 
major components: 

P(T) = Po + Pph(T) + pso(T), (Bl) 
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n 

1 2 

3 4 

</n (10 4 dyn/cm) 

40.905 7.402 

-1.643 -0.609 

<t>n (10 4 dyn/cm) 

16.685 -4.051 

3.267 0.424 

4>n (10 4 dyn/cm) 

9.616 -0.841 

0.603 -0.501 

el( u (10 4 dyn/cm) 

0.000 0.632 

0.000 -1.092 

tn l (10 4 dyn/cm) 

0.000 -0.632 

0.000 -1.092 


TABLE I. Parameters of the 4NN force constant model with 
the off-diagonal couplings. The corresponding dispersion re¬ 
lation is shown in Fig. [6] 



FIG. 7. Comparison between our simulations (green dotted 
lines) and experimental data? 7 (red markers) for the resis¬ 
tivity of graphene at carrier density varying from 1 (top) to 
3 x 10 12 cm -2 (bottom) with a step of 0.5 x 10 12 cm' 2 . The 
blue horizontal dotted lines represent the impurity-limited re¬ 
sistivity po- 



HC. 8. p p h (blue makers) and pso (red markers) versus 
temperature at electron densities 1 (empty markers) and 
13.6 x 10 12 cm' 2 (filled markers). 


Appendix C: Band structures of CNTs 

Figure [9] shows the band structure of metallic armchair, 
metallic zigzag and semiconducting zigzag CNTs for se¬ 
lected diameters. For small diameter CNTs, only the 
lowest band is involved in the transport. The effective 
mass of semiconducting CNTs decreases with increasing 
diameter, leading to higher mobility in CNTs with large 
diameter. 

Figure [TO] shows the position of the Fermi level in the 
band structure of CNTs with diameter of 3.5 nm. This 
figure helps to interpret the results presented in Fig. [2]). 
At high carrier density, more bands are included in the 
transport energy window. That not only brings in bands 
with finite effective mass, but also more scattering mech¬ 
anisms. 


where T is temperature, po is due to impurities, p p h is 
due to the intrinsic phonons of graphene, and pso is due 
to the surface optical phonons of the substrate. p p h is 
calculated as described in the main document, pso is 
calculated with the model of Ref. in which Si0 2 is 
considered for the substrate, po is determined by taking 
the difference between experimental data and the sum of 
p p h and pso at 30 K. Figure [7] compares the results of 
the calculation with experimental data^ 7 in a wide range 
of temperature and for different electron densities n 2 p. 
It is clear that they agree very well with each other. 

Figure [8] shows the temperature dependence of p p h 
and pso- At l° w temperature, the effect of the surface 
optical phonon scattering is not significant because the 
lowest surface optical phonon energy considered in the 
model is 59.98 meV. The surface optical phonon scatter¬ 
ing is important at high temperature when the carrier 
density is low. It almost dominates the resistivity at 
ri 2 D = 10 12 cm' 2 . However, at higher densities, its con¬ 
tribution becomes negligible due to strong screening.— 


Appendix D: Carrier density dependence of p x U 2 d 

Figure [Tl] shows that, at 300 K, the product p x n 2 D 
does not depend on the carrier density n 2 £> in graphene. 
The same behavior is found for CNTs at high carrier den¬ 
sity but only approximately, there are deviations coming 
from band structure effects. 
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FIG. 9. Band structures of metallic armchair (a-d), metallic zigzag (e-h) and semiconducting zigzag (i-1) CNTs with diameter 
of 3.5 (a,e,i), 6 (b,f,j), 10 (c,g,k) and 16 nm (d,h,l). The 2D effective carrier density is fixed at 10 12 cm -2 . The horizontal line 
represents the Fermi level at 300 K. 
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FIG. 10. Evolution of the Fermi level at 300 K in the band structure of armchair (a-d), metallic zigzag (e-h) and semiconducting 
zigzag (i-1) with diameter of 3.5 nm and 2D effective carrier density of 1 (a,e,i), 6 (b,f,j), 13 (c,g,k) and 20 x 10 12 cm -2 (d,h,l). 
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T — 300K, d ss 3.5 nm 
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FIG. 11. Mobility (fj,) multiplied by the carrier density (tl 2 d) 
versus ri 2 D in graphene and CNTs (T = 300 K, d « 3.5 nm). 
The same data but presented as p versus U 2 D are shown in 

Fig.gb. 
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